Massive gravitons trapped inside a hypermonopole 



00 ■ 

o : 
o . 

(N ■ 

m ; 



i ' 

<d : 

43 . 



> 

q 

ON 

o 

00 

o 



X 



Antonio De Felic^B and Christophe Ringevalll 
Theoretical and Mathematical Physics Group, Centre for Particle Physics and Phenomenology, 
Louvain University, 2 Chemin du Cyclotron, 1348 Louvain-la-Neuve (Belgium) 
(Dated: September 2, 2008) 

We propose a regular classical field theory realisation of the Dvali-Gabadadze-Porrati mechanism 
by considering our universe to be the four-dimensional core of a seven dimensional 't Hooft-Polyakov 
hypermonopole. We show the existence of metastable gravitons trapped in the core. Their mass 
spectrum is discrete, positive definite, and computed for various values of the field coupling constants: 
the resulting Newton gravity law is seven-dimensional at small and large distances but can be 
made four-dimensional on intermediate length scales. There is no need of a cosmological constant 
in the bulk, the spacetime is asymptotically flat and of infinite volume in the extra-dimensions. 
Confinement is achieved through the local positive curvature of the extra-dimensions induced by 
the monopole-forming fields and for natural values of the coupling constants of order unity. 

PACS numbers: 04.40.-b, 04.50.-h, ll.10.Kk, 98.80.Cq 



INTRODUCTION 

Gravity occupies a central role in high energy physics 
and cosmology. On one hand, the unification of the 
fundamental interactions in the context of String The- 
ory suggests that we may live in a more than four- 
dimensional world [j], Q. On the other hand, the re- 
cent acceleration of our universe has been confirmed by 
different experiments and it is now a widely accepted 
important result of modern observational cosmology 
The idea that such an unexplained acceleration may be 
the signature of extra-dimensions has been intensively 
explored in the recent years 0, [H, @j- In the Dvali- 
Gabadadze-Porrati (DGP) model, the extra-dimensions 
(bulk) may actually be non-compact and of infinite vol- 
ume [3, [i| . Gravitons are reflected back onto our universe 
(brane) due to a different gravity coupling constant on 
the brane and in the bulk. The original DGP action in 
n c + 4 dimensions reads 



\Vg\Rd 4 x 



g\Rd n ^X, (1) 



where g and R are respectively the determinant and 
scalar curvature of the induced metric along our brane, 
while g and R are the corresponding quantities in the 
bulk. It has been shown that this model could actu- 
ally explain the observed acceleration of the universe, 
although some works suggest that it may be spoiled by 
instabilities 

@, EU E3- Although the form of Eq. © 
has been originally explained by quantum effects, "reg- 
ularised models" have been proposed to justify it from 
a more classical and tractable point of view, free of 
instabilities. Such an approach has been explored in 
Refs [H, EH, EH an d shown to confine gravitons by ex- 
plicitly choosing some profile for M*(X) or g(X). In 
a complete physical framework, both of these functions 
are however not free and it is not clear that the DGP 
mechanism could indeed appear in any classical system. 



This question is of crucial importance in order to assess 
the viability of both infinite volume extra-dimensions and 
instability-free DGP-like mechanism. 

In this letter, we answer this question in the con- 
text of canonical classical field theory. Our approach is 
motivated by condensed matter physics: topological de- 
fects are a direct consequence of the symmetry breaking 



mechanism and can model smooth branes 15|, llq . |17 |. 
Assuming the space-time to be seven-dimensional, an 
SO (3) spontaneous symmetry breaking in n c — 3 codi- 
mensions generically forms 't Hooft-Polyakov hyper- 
monopoles In the following, we prove the ex- 

istence of a DGP-like mechanism in the core (assumed to 
be our universe) of such a monopole. 

Compared to lower dimensional defects [201 ]. the ex- 
istence of positively curved n c — 1 dimensional regions 
in the bulk is crucial to allow metastable gravitons to 
be trapped inside the core. Six is indeed the minimal 
number of spatial dimensions for which there exists a fo- 
liation of the extra-dimensions by two-dimensional pos- 
itively curved surfaces. In order to allow for a varying 
Planck mass, we have for completeness included a dila- 
ton ip having a mass ma in the Einstein frame. In the 
Jordan frame, the action associated with this system is 

1 
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AB 



(2) 



where the dilaton potential reads U = m^ip 2 exp(2?/>/5). 
The 50(3) Higgs field 3? = {4> a } is in the triplet repre- 
sentation (a S {1,2,3}). Its vacuum expectation value 
v breaks 50(3) into U(l). The covariant derivatives T>a 
enforce gauge invariance and incorporate the gauge fields 
Ca = {C%}, 

V A $> = d A ® - qC A A (3) 
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where r, 8, <p are spherical coordinates in the three extra- 



FIG. 1: Field and metric profiles forming the hypermonopole 
for a = 2.05, e = 0.50 and = 1.00 (top). The space-time 
is flat asymptotically and in the core, but strongly curved in 
the intermediate region. The spatial sections for 9 — n/2 are 
represented in the bottom panel as a function of the radial 
coordinate. 



q being the charge, while the field strength tensor H ab 



H A b = d A C B -d B C A - qC A A Ci 



(4) 



As for the dimensional analysis, we have [k 2 ] = M 5 , 
[q] = M- 3 / 2 , [A] = M- 3 and [C A ] = [$] = [v] = M 5 / 2 . 



BACKGROUND GEOMETRY 

Static self-gravitating monopole configurations associ- 
ated with the action @ can be obtained by imposing 
isotropy in the extra-dimensions, plus Poincare invari- 
ance along the four internal brane coordinates x^. Our 
ansatz for the metric is 



dimensions and dfr 



d9 2 



For the 



monopole-forming fields, the internal space of SO(3) is 
mapped to the three extra-dimensions with a purely ra- 
dial Higgs field 



* = vf(r)u r 



and winding gauge fields 



1 - Q(r r 
C e = u 

q 



1 - Q(r) . 
C„ = smoue 



(6) 



(7) 



the other components vanishing. Here, f(r) and Q(r) are 
two dimensionless functions such that, far from the core, 
f(r) — > 1 and Q(r) — > to recover a Dirac monopole. 
In the core, regularity imposes /(0) = and Q(0) = 1. 
Concerning the metric coefficients, the energy associated 
with the defect being finite and localised, we look for 
asymptotically flat spacetime, a — > 0, oj — > r and ip — ► 0. 
Regularity in the core also imposes cr'(0) = ip'(0) = 
and u/ ~ r. 

The system of coupled non-linear differential equations 
obtained from the action is of order ten and does not 
have any obvious analytical solution. Once the radial 
coordinate is expressed in unit of the Higgs Compton 
wavelength, the differential system is parametrised by 
three dimensionless parameters 



2 2 

a = k v , 



2 2 

q v 



Xv 2 



(8) 



d^ 



^ri^dx^dx" + dr 2 + w(r) 2 dfi 2 , 



(5) 



where and are respectively the mass of the Higgs 
and gauge bosons. Under the above-mentioned bound- 
ary conditions, the numerical integration of the equations 
of motion is a challenging problem that has been over- 
come by using recent advances in the field 2l|. We have 
found monopole solutions for almost any values of the 
above parameters; only when the stress energy becomes 
super-Planckian the system develops some singularities 
preventing the spacetime to be asymptotically flat. As 
can be seen in Fig. [U the Higgs and gauge field profiles 
are typical of topological defect configurations while the 
dilaton is gravitationally trapped inside the core. The 
profile of a(r) traces the gravitational redshift: clocks 
are ticking differently inside and outside the monopole. 
More interesting is the profile of uj(r). Up to a 4w fac- 
tor, u> 2 (r) gives the area of the two-sphere of radius r 
in the extra-dimensions. As can be seen in Fig. [TJ there 
is a region at finite distance from the core where u>(r) 
does no longer grow as r but remains almost station- 
ary: the extra-dimensions become cylindrically shaped. 
As we show in the next section, gravitons become reso- 
nant at these length scales and metastable from a four- 
dimensional point of view. Notice that the spacetime is 
non-compact and asymptotically Minkowski. 
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TENSOR FLUCTUATIONS 

We now consider the four-dimensional tensor perturba- 
tions around the previously computed background. The 
perturbed metric is given by Eq. (|5|) upon the replace- 
ment ?7p„ — > + hf^u, where is a spacetime de- 
pendent transverse and traceless tensor. The linearised 
equations of motion for are obtained by expanding 
Eq. d2|) at second order and have already been derived 
for an arbitrary number of extra-dimensions in Ref . [20| . 
Defining the dimensionless conformal radius z(r) and ten- 
sor £ 



as 



m h 



exp(-cr/2)dr, £ 



e V-/2 e 3a/4 wft 



(9) 



the equation of motion for the spin-two fluctuations can 
be recast into 



d 2 £ 



w 2 + w 



IS - □ £ = 0, (10) 



where the tensor indices have been omitted. Deriva- 
tives are with respect to z, □ = m^ 2 r] fJ -' y d t _ l d l/ is the 
d'Alembertian along the brane, and 



W 



UJ 



do 



d 2 



tan I 



fill 



After a four-dimensional Fourier transform on the brane 
coordinates, and an expansion over the spherical harmon- 
ics in the angular extra-dimensions, we have in unit of the 
Higgs mass □ -> -Tj^p^p" = M 2 and L 2 -> -£(£ + 1). 
One immediately recognises in Eq. (|10p the Schrodinger 
equation of a supersymmetric quantum mechanical sys- 
tem in a central potential V 2 = W 2 + W [22]. The 
operator L 2 is the angular momentum, W(z) is the su- 
perpotential and M 2 plays the role of the energy. A 
subtlety is that our coordinate z lies on the positive axis 
only. However, since h^v must remain finite on the brane, 
Eq. (JHJ) implies that £ M „ should vanish in z — 0. Under 
this condition, and the usual normalisability at infinity, 
the differential operators remain regular enough to use 
the results of supersymmetric quantum mechanics. For 
M 2 = L 2 = 0, Eq. (JTUJ) is solved by the "ground state" 
£o oc wexp(3cr/4 -1- ip/2) which is however not normalis- 
able asymptotically. The ground state of the superpart- 
ner potential V\ = W 2 — W is l/£o, which is not regular 
in z — 0. As a result, "supersymmetry" is broken and the 
spectrum is necessarily positive definite, M 2 > 0: there 
is no massless mode neither tachyon on the brane. 

Introducing the orthonormal basis of eigenfunctions 
UM,e{z), solutions of Eq. (fT0|) , such that 

u* M4 (zi)u M A z 2)dM = 5(zi - z 2 ), (12) 

one can check that the retarded Green function for £ can 



— v, 


= w 2 


+ w 


— v , 


= w 3 - 


w 




FIG. 2: Superpartner potentials and Vi(z) confining 

the metastable gravitons. The background fields are those of 
Fig. [I] 



be expanded as [23| 
Gt(X Vt X 2 ) = - j ^e^K-^)^ rr( 1;V?1 ) 

um,i{zi)u* m e (z 2 )dM 



(27T) 



M 2 + (p) 2 - (p° + ief ' 
(13) 

Using the above equation together with Eq. , the ten- 
sor modes sourced by any transverse and traceless stress 
tensor S^ U {X) are given by 



M*0 



2k 2 



_ e -V(2l)/2 e -3a( Zl )/4 



x / G 5 (A i; A 2 )e-^( Z2 )/ 2 e +3<T ( Z2 )/ 4 W (z 2 )5^(A 2 )d 7 A 2 . 

(14) 

In order to gain some intuition on the previous expres- 
sions, let us first consider the case of a seven-dimensional 
flat spacetime. Setting ip = a = everywhere, as well 
as ui — r, Eq. (I10| can be integrated and the normalised 
modes are 



' t+1/2 



(Mz). 



(15) 



Considering a "point-like" static source s fiU (x) on the 
brane 



Sfj,u(X) 



lim — 

z^O z 



: S(z)S (cos 0) 6(tp) Spu (x) , 



(16) 



we can explicitly integrate Eqs. (fTB")) and (|14[) using the 
flat modes u M e . Only the s-waves {I = 0) have a non- 
vanishing contribution on the brane since for z — > 0, one 
has u b AI £-£ (z) J z — > 0. The four-dimensional integral over 
the momentum p in Eq. (|13j) , together with the denom- 
inator containing M 2 , is the classical retarded Yukawa 
propagator. After some calculations, one finally gets on 
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FIG. 3: Spectral density p(M) as a function of the graviton 
mass M. There are at least two trapped metastable gravitons 
in the background geometry of Fig. [T] The long lived reso- 
nance is well fitted by a Breit-Wigner distribution centered 
at m g ~ 0.10326 with a width r ~ 3.5 x 10~ 4 (in mh units). 



the brane 



hu. v (xi) = lim 



K f 
87T 2 TOJ J 



'M. 



(17) 



I Ail 



where Ax = x\ — x 2 . Using the expansion of the Bessel 
'2/nMz, the previous expression 



function, 
simplifies to 



//// 



47T 3 7 



X 2 



J fXV 



(X2) 



\Ax\ 



(18) 



which is the standard linearised solution of the Einstein 
equations around a seven-dimensional Minkowski space- 
time. Notice the power law dependence l/|Aa?| d ~ 2 in 
d = 6 spatial dimensions, as well as the 4-7T 3 factor which 
is d — 2 times the surface of the unit (d — l)-sphere, as 
one would have obtained from the Gauss law. 

In the background geometry of the hypermonopole, 
the situation is nearly the same apart that the mode 
functions umj{z) are now modified. As can be seen 
in Eq. (fT7|) . we need the values of the rescaled spectral 
density associated with s-waves on the brane p(M) — 
\um o(0)| 2 /I u m o(0)I 2 - The tensor modes are then given 
by ' 



2 K 2 e-^°> 
87r 3 m 2 



£{p(M)M 2 } 
d x 2 — - — — — - 



s^(x 2 ), (19) 



where £{■} stands for the forward Laplace transform 
evaluated at |Aa?|. From the monopole- forming fields 
computed in the previous section, we have plotted in 
Fig. [2] the potential V 2 (z) and its superpartner. Since 
the spacetime is asymptotically flat, V 2 vanishes at in- 
finity and there is not any bound state. However, V 2 



(and also V\) exhibits a barrier at the location of max- 
imum curvature allowing metastable modes in the core. 
We have numerically solved the equation of motion (Ti"0|) 
for these potentials and plotted in Fig. [3] the resulting 
spectral density. When this quantity is constant, grav- 
ity is purely seven-dimensional on the brane. This is the 
case for large, but also for low values of M 2 , as in the 
Gregory-Rubakov-Sibiryakov model [24| . 

Notice that the constant value of p gives the effective 
gravitational coupling constant: here, it is different for 
large and low values of M due to the dilaton condensation 
in the core as well as the gravitational redshift. In the in- 
termediate range, p(M) is strongly peaked for particular 
values of M: these are the resonant metastable modes. 
We have also numerically checked that there is not any 
bound state with M 2 < 0, as expected from the super- 
symmetry arguments. Changing the background param- 
eters a, (3 and e affects the mass spectrum and the width 
or number of metastable modes can be easily adjusted. 
Lowering e delocalises the gauge fields (in unit of the 
Higgs Compton wavelength) and the position of the bar- 
rier is pushed towards larger values of z. Decreasing a, 
or increasing e, reduces the height of the barrier while (3 
changes its shape. Let us notice that if a is too small, 
or e too big, we do no longer observe resonances and this 
corresponds to the disappearance of the confining nest 
on V 2 . This is reminiscent with the properties of bound 
states in the case of broken supersymmctry. 

To understand how these resonances realise the DGP 
mechanism, one can approximate them as Dirac distri- 
butions. Let us say we have one trapped graviton at a 
mass m g , then neglecting the smooth changes in the spec- 
tral density, p(M) ~ 1 + C5(M — m g ) where C encodes 
how peaked the resonance is. The Laplace transform in 
Eq. (fl9|) simplifies to 



£{M 2 p(M)} 



I A: 



13 + Cm e e ~ 



(20) 



and four-dimensional gravity is recovered over the length 
scales (rng/C) 1 / 3 < |Ax|m g < 1, provided the mode is 
light enough m g < C. At small and large distances, 
seven dimensional gravity is recovered while in between 
we have even observed some fractional power dependen- 
cies. When more than one gravitons are trapped the 
situation becomes even more complex and the detailed 
analysis of these effects is left for a forthcoming work. 



CONCLUSION 

We proposed here a canonical classical field theory 
model which describes a seven-dimensional monopole, at 
the core of which gravitons get trapped. Their mass spec- 
trum being positive definite, there are no instabilities for 
the tensor modes. This phenomenon turns out to be a 
natural way, in the context of field theory, to implement 
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the DGP idea. The required field configurations can be 
obtained without fine-tuning from a dense set of coupling 
constant values of order unity. However, to obtain a four- 
dimensional gravity behaviour over a wide range of length 
scales, some amount of fine-tuning is certainly required 
to confine an almost massless mode. Notice that the 
smaller the graviton mass, the larger the effective four- 
dimensional Planck mass, possibly addressing the mass 
hierarchy problem. It would be interesting, if possible, to 
find a condensed matter system for which this mechanism 
could be experimentally explored. 
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